Realistic quantum manipulation of two-level system fluctuators 
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Two-level system fluctuators in superconducting devices have demonstrated coherent coupling 
with superconducting qubits. Here, we show that universal quantum logic gates can be realized in 
these two-level systems solely by tuning a superconducting resonator in which they are imbedded. 
Because of the large energy separation between the fluctuators, conventional gate schemes in the 
cavity QED approach that are widely used for solid-state qubits cannot be directly applied to the 
fluctuators. We study a scheme to perform the gate operations by exploiting the controllability of the 
superconducting resonator with realistic parameters. Numerical simulation that takes into account 
the decay of the resonator mode shows that the quantum logic gates can be realized with high fidelity 
at moderate resonator decay rate. The quantum logic gates can also be realized between fluctuators 
inside different Josephson junctions that are connected by a superconducting loop. Our scheme 
can be applied to explore the coupling between two-level system fluctuators and superconducting 
resonators as well as the coherent properties of the fluctuators. 
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I, INTRODUCTION 

Spurious two-level system (TLS) fluctuators are con- 
sidered a serious source of low-frequency noise in su- 
perconducting qubitsA, and the characterization of these 
fluctuators in solid-state devices has a long history 2 . 
Most recently, coherent coupling between TLS fluctu- 
ators and a superconducting phase qubit was observed 
via the novel energy splittings in spectroscopic measure- 
ments^^. It was shown that the TLS fluctuators have 
much longer decoherence times than the superconduct- 
ing qubits, raising the possibility of realizing quantum 
manipulation on these fluctuators^. 

The key question in manipulating the TLS fluctua- 
tors is how to implement the required coherent manip- 
ulation and readout. Located sparsely inside solid-state 
devices, the fluctuators usually do not interact with each 
other, and their states are hard to control. The coupling 
between the fluctuators and solid-state devices provides 
us with a tool to achieve the quantum manipulation^. 
However, conventional gate schems using cavity QED ap- 
proach that are usually exploited for solid-state qubits 
cannot be applied to this system because of the large en- 
ergy separation between the fluctuators. In this work, we 
will present a gate scheme that exploits the controllabil- 
ity of the superconducting resonator to implement high 
fidelity gates on the TLS fluctuators ^ 8 ^ 10 ' 11 , even when 
the decay of the resonator is a few megahertz. The super- 
conducting resonator acts as a knob that controls the dy- 
namics of individual fluctuator, as well as coupling them 
together. Working with practical parameters from the 
superconducting Josephson junction resonator, we will 
design single-qubit and two-qubit quantum logic gates 
in the presence of resonator decay. Our scheme takes 
into account the full coupling Hamiltonian between the 
TLS fluctuators and the resonator. Readout of the fluc- 
tuators can also be performed by measuring the trans- 
mission through the resonator. This scheme can be ex- 
tended to fluctuators in different Josephson junctions by 



connecting the junctions into the same superconducting 
loop due to the nonlocal nature of the microwave mode 
of the resonator. This work hence provides a realizable 
design for coherent manipulation of multiple TLS fluc- 
tuators, which is closely related to current experimental 
efforts in studying the fluctuators and their coupling with 
superconducting resonator modes. 

Various superconducting resonators in the microwave 
regime, including superconducting transmission lines, 
Josephson junctions, SQUID's, and superconducting 
lumped element resonators, have recently been demon- 
strated and have shown quantum behavior and strong 
coupling with superconducting qubit a 12 ' 13 ! 14 ' 15 . Super- 
conducting resonators are also promising systems for 
studying quantum effects such as single photon gener- 
ation and lasing— and one of us has shown recently that 
a Josephson junction can be used to probe various prop- 
erties of TLS fluctuators, e.g. to resolve the mechanism 
that couples the fluctuators to the junction^. While we 
will focus on the Josephson junction resonator, we want 
to emphasize that our results can be readily generalized 
to other superconducting resonatorsA^&. The paper is 
organized as the following. In Sec. HH we will study the 
coupled system of the fluctuators and a Josephson junc- 
tion resonator, including the driving on the resonator. 
In Sec. IIII| we will derive the effective Hamiltonian for 
the TLS fluctuators in the dispersive regime where the 
quantum operations are implemented. We will also de- 
rive the residual coupling between the fluctuators and 
the resonator in this regime. In Sec. \W\ we will present 
detailed scheme for single-qubit and two-qubit quantum 
logic gates. Then, we will estimate the decoherence of 
the fluctuators during the gate operations in Sec. [Vl We 
will also test the fidelity of the quantum operations with 
numerical simulation of the full Hamiltonian, taking the 
resonator decay into account. In Sec. Wl\ we will dis- 
cuss the readout of the fluctuators and the extension of 
gate scheme to fluctuators inside different junctions. The 
conclusions will be given in Sec. I VIII 
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Figure 1: A Josephson junction resonator containing spuri- 
ous two-level system fluctuators denoted by arrows, (a) Fluc- 
tuators in a single junction and (b) fluctuators in different 
junctions. 



II. THE SYSTEM 

Consider the system in Fig.[T](a), where TLS fluctua- 
tors inside the amorphous layer of a Josephson junction 
couple with the junction resonator in an RF SQUID loop. 
With total capacitance Co, Josephson energy Ej, loop 
inductance L, and magnetic flux <& ex inside the SQUID 
loop, the Hamiltonian of the resonator can be written as 



p2 
2C7n 



- Ej cos(2e$/fr) 



2L 



(1) 



in terms of the phase <I> and the conjugate momentum P$. 
This Hamiltonian can be approximated as an oscillator 
mode with a phase shift & s from the origin: 



H c « Pj/(2C ) + C> c 2 ($ - $ s ) 2 /2, 
and the phase shift satisfies 

H$ s + 2eLE j sin(2e$ S /H) = -H$ ex . 
The frequency of the resonator can be written as 



Ae 2 Ej cos(2e<f> s /h) 



LC 



(2) 



(3) 



(4) 



which can be tuned in a large range by the magnetic flux 
&ex- In addition, driving can be applied to the resonator 
by e.g. applying an external radio-frequency current SI C 
to the resonator with 5I C &. 

The TLS fluctuators reside inside the tunneling 
layer and can couple with the junction resonator by 
various mechanisms. For example, the coupling to 
the critical current of the junction takes the form 
— (2e/h)EjQ J2 n in ■ <?„, where j n is the polarization and 
magnitude of the coupling. Denoting t he resonato r anni- 
hilation operator by a with $— <£> s = y/h/ (2Cow c )(a+a T ). 
Let j n = (j xn , 0,0) for simplicity and ujd be the driving 
frequency. The total Hamiltonian of the coupled system 



in the rotating frame can be written as 
Hi 

Hi 



H c + H x + H K 
hAcd^a + e(a + a') 



(5) 
(6) 



y"l [{hA n /2)a nz + g n (aa n+ + aV„_)] (7) 



H K = 2J hajkajak + c k {a\a + a^a k ) 



(8) 



where H c is the Hamiltonian of the driven resonator 
mode with the det uning A c = lu c — LUd and the driving 
amplitude e = SI c ^/h/(2CouJ c ), Hi is the Hamiltonian of 
the fluctuators including the coupling between the fluctu- 
ators and the resonator mode, and H K is the Hamiltonian 
of the thermal bath connected to the resonator. Here, the 
index n labels different fluctuators, cr na are the Pauli op- 
erators, A„ = Lo n — LUd is the detuning of the fluctuators, 
and 



9n = E J j xn y / h/(2C uJc) sin(2e$ s /ft) 



(9) 



is the coupling constant. Note that coupling constant 
for other coupling mechanisms such as dielectric coupling 
between the fluctuators and the resonator can be derived 
similarly^. The decay of the resonator is modeled by its 
coupling to a bath of modes described by the annihilation 
operator a k with frequency u>k and coupling constants Ck ■ 
The decay rate is given by k = '^J2 c k^( LJ ~ ^fc}^- The 
Hamiltonian H t describes a typical cavity QED system 
between the fluctuators and the junction resonator—. 

Note that the driving on the resonator generates a 
time-dependent oscillation in the phase variable with the 
amplitude 6§d — 5I c /Cqlo 2 . To keep the nonlinear term 
in the Josephson energy to be small, the oscillation ampli- 
tude needs to be small, e.g. |2e£$d/fi| < 0.1. With 1/L ~ 
Ae 2 Ej/h 2 and typical parameters Ej ~ 27rx 100 GHz and 
Co ~ 10 _12 pF, we estimate that the driving amplitude 
is bounded by e < 27r x 1GHz. 



III. THE DISPERSIVE REGIME 

In this work, we study the quantum logic operations 
in the dispersive regime where the coupling g n is much 
weaker than the detuning between the fluctuators and 
the resonator: g n <C |A nc | with A„ c = A„ — A c . In this 
regime, we can apply the following unitary transforma- 
tioni 



U 



= e - £ (o-ot)/A c 



n- 



-S„(at 



h a)/A„ 



(10) 



to the system. After the transformation, the Hamiltonian 
becomes H t = UH t U^ with H t = H c + Hi + H x to the 
second order of g n /A nc . The Hamiltonian is now divided 
into three parts: a Hamiltonian for the resonator H c , 
an effective Hamiltonian for the fluctuators Hi, and a 
small residual coupling between the fluctuators and the 
resonator H x . 
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The Hamiltonian Hi can be written as 



H, = 



E 



a. 



H mt + H k (ll) 



2 2 

Hint = y \nn( (J n+< J m- + <Tm+ a 'n- 

)/2 (12) 

H K = ^ j (g n c k /& nc ){(j n+ a k + a { k (j n -) (13) 

n, k 

which includes the effective single qubit terms, an 
exchange- like interaction H int , and an induced coupling 
to the bath modes of the resonator H K . We derive the 
detuning for the single qubits as 

A„ = A„ + ( 5r 2 i /A„ c )(l-2e/A c ) (14) 

and the Rabi frequency as 



(15) 



The coupling constant in the exchange-like interaction 
can be derived as 



^mn — ym^7n(A mc -(- A nc )/ (A mc A nc ; 



(16) 



In the following section, we will study the implementation 
of the quantum logic gates with the Hamiltonian Hi . 
The residual coupling H x can be written as 



H, 



a) a + e 



A c - 2A f , 

2 A nc A c 



(a + a f ) 



(17) 

where the first term is the Stark shift for the resonator 
and the second term is a coupling to the resonator am- 
plitude originated from the finite driving amplitude. Be- 
cause of the amplitude shift in the unitary transforma- 
tion in Eq. (flO|) , the average occupation of the resonator 
is now zero with (a^a) w 0. Hence, the first term has a 
small effect on the fluctuators during the quantum oper- 
ations. The second term can induce a small modification 
to the coupling constant X mn in the effective interaction 
in Eq. (fl2|l which will be studied in detail below. 



IV. QUANTUM LOGIC GATES 

Universal quantum gates can be performed by control- 
ling the effective Hamiltonian Hi. Here, we present the 
scheme for the single-qubit and two-qubit gates with typi- 
cal parameters from superconducting Josephson junction 
resonators. 

Single-qubit gates on a chosen TLS (e.g. n = 1) can 
be performed by adjusting the frequency of the driving 
source to be close to the frequency of this TLS. By ad- 
justing the driving amplitude and the detuning of the 
resonator, the effective qubit parameters A x and Qi x can 
be adjusted in a wide range. For a given detuning A c , 
one can adjust the driving amplitude e in Eq. lfl4|) to 




200 
A c (MHz) 

Figure 2: Upper plot: gate times for Hadamard gate (solid 
curve) and spin flip gate (dashed curve) versus detuning A c . 
Lower plot: effective couplings (32 and (3' 2 versus A c (main 
plot); energy E\ (solid curve) and -E2 (dashed curve) versus e 
at A c = 300 MHz (inset). 





A c (2ttxMHz) 


e (2ttxMHz) 


Qi x (2ttxMHz) 


time (ns.) 


X 


120 


-60 


60 


8.3 


H 


160 


-32 


21.3 


16.6 



Table I: Example parameters for implementing the spin flip 
gate (X) and the Hadamard gate (H). 



have Ai = and obtain the spin flip gate X. The gate 
time can be found to be r g = 7rA„ c /2eg„ at the chosen 
e. One can also adjust e to have Ai = fli x to implement 
the Hadamard gate H. Here, the driving amplitude is 
chosen to be 



(A e A„ e + gl)A c 
2(A C + g n )g„ 



(18) 



with the gate time r g — irA nc /2V2eg n . In Fig. [21 
we plot the gate times of the Hadamard gate and spin 
flip gate with the parameters Ai = 2ir x 40 MHz and 
gi = 2tt x 40 MHz. It can be shown that gate times on the 
order of 10 ns. can be achieved. In Tabled we list two sets 
of gate parameters at the detunings A c = 2-7T x 160 MHz 
for the Hadamard gate and A c = 2ir x 120 MHz for the 
spin flip gate respectively as an example. The corre- 
sponding driving amplitudes are listed in the table. For 
fluctuators not involved in the gate (e.g. n = 2), we have 
hA n ^> Q nx due to the fact that the fluctuators are well 
separated in energy. Hence, the gate operation only in- 
duces a dynamic phase to these fluctuators. Meanwhile, 
the effective interaction between the fluctuators is also 
prevented from generating controlledj)hases due to the 
off-resonance condition Ai„ <C |Ai — A„|. 
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Two-qubit gates can be performed via the effective 
exchange- like coupling in Eq. (fTTj) . This coupling can 
generate SWAP gate and V SWAP-like controlled gates 
when the two qubits are near resonance with A m — A„ ~ 
022. However, as noted above, the TLS fluctuators in 
our system are usually far off-resonance from each other 
due to the large energy separations and are also hard 
to be manipulated individually. We will now show that 
two fluctuators can be brought into effective resonance 
by controlling the resonator mode, and hence high fi- 
delity two-qubit gates can be performed. We first rewrite 
the single TLS energy in Hi as ^ n (£„/2)ff„ z where 



El = A? 



QJL and 



(19) 



with co&9 n — A n /E n and sin( 



tt nx /E n . Note <r„ 



and a ny can be defined similarly. As the driving ampli- 
tude increases, the Rabi frequency increases accordingly 
and the detuning A„ will be affected. As plotted in the 
inset of Fig. [21 a driving amplitude can be found where 
E n is the same for the two fluctuators. Denoting these 
fluctuators as n = 1 and n — 2, we have E\ = E2 at this 
point, and they are now in resonance. In the interaction 
frame, the effective Hamiltonian for these two fluctuators 
then becomes 

H12 = 0\Oi z a 2z + 2 {a x+ d 2 ^ + ai^a 2+ ) (20) 

with the coefficients 



02 = (A 12 /4)(1 + A 1 A 2 /S 1 2 ) 



(21) 
(22) 



respectively. Meanwhile, at large driving amplitude with 
e ~ A c , the second term in the residual coupling H x 
in Eq. lfl~7|) induces virtual transitions between the res- 
onator and the fluctuators, which modifies the coupling 
constant 2 to become 0' 2 . Denoting the second term in 



Eq.|l7} as f n a T , 



la 



02=02 + 



r), we can derive 

A/ 2 (£i + E 2 - 2A C ) 
2(E 1 -A C )(E 2 -A C )' 



(23) 



In Fig. [2 we plot both the couplings 2 and 0' 2 versus 
the resonator detuning A c for comparison. A small but 
finite modification of the coupling coefficient can be seen. 
Note when \E\ y2 — A c | ~ Ai 2 , the second term in H x 
induces real transitions between the fluctuators and the 
resonator, which can seriously affect the gate operations 
and should be avoided when designing the gates. 

Two-qubit gates of the form of So = exp(— ii/[ 2 l t) can 
now be performed in the rotating frame. A SWAP gate 
has the gate time r g = tt/2\0' 2 \. The 0\ term in Eq. lj20|) 
contributes only phase factors in the computational basis 
to the gate operation. With the following parameters: 
Ax = 0, A 2 = -2tt x 60 MHz, 9l = 2ir x 40 MHz, g 2 = 
2tt x 30 MHz, and A c = 2tt x 300 MHz, we find that E 1 = 





1-Qubit 


2-Qubit 


swap op. 




ivA nc /g n € 


tt/2|^| 


n/2g n 


-1 
T d 


2 /a 2 


2 /a 2 


k/2 


Tg/Td 


10~ s 


0.01 


0.02 



Table II: Gate times r 9 , decoherence rate r~ , and the ratios 
Tg/T d for single-qubit and two-qubit gates. The column la- 
belled as "2-Qubit" is for the two-qubit gates studied in our 
scheme and the column labelled as "swap op." is for the swap 
operation in the Cirac-Zoller gate^. 



E 2 = 2tt x 74.0 MHz at e = 2tt x 277.2 MHz. Here, 0' 2 = 
—2tt x 1.8 MHz and the SWAP gate can be performed 
with a gate time of r g = 137.9 ns.. 

We note that controlled quantum logic gates can also 
be performed using the Cirac-Zoller gate which was first 
studied in ion trap quantum computing^. This gate in- 
cludes three steps: a swap gate between the first TLS 
and the resonator, a conditional phase gate between the 
resonator and the second TLS, and another swap gate 
between the first TLS and the resonator. The swap gate 
in the first and third steps can be implemented by choos- 
ing A c = Ai and e = 0, and the gate time is 7r/2<7i 
which is of the order of a few nanoseconds. In the con- 
ditional phase gate, the driving frequency is close to the 
chosen qubit but is still in the dispersive regime. The 
Stark shift (g 2 / 'A 2c )a 2z a^ a generates a conditioned phase 
shift on this TLS when the resonator is in state |1). The 
gate time is t cg = nA 2c /2g 2 . At g 2 — 2tt x 30 MHz 



and A 2c = 2tt x 120 MHz, t cg - 30 ns.. Note that the 
first TLS is subject to stronger decoherence during the 
swap operation due to its near-resonance coupling with 
the resonator, as will be discussed below. 



V. DECOHERENCE 

The intrinsic decoherence of the TLS fluctuators is very 
slow and can be ignored during the gate operation. How- 
ever, the coupling between the resonator and the fluctua- 
tors induces decoherence that cannot be neglected. In the 
dispersive regime, the decoherence rate can be calculated 
from the noise term H K in Eq. lfTTT) . It can be shown that 
the decoherence rate is on the order of t^ 1 ~ g^nj A 2 nc 
during the quantum logic gates with t^ 1 -C k in the dis- 
persive regime. In contrast, the decoherence rate during 
the swap operation in the Cirac-Zoller gate is t^ 1 ~ «/2 
which is much faster than the decoherence rate in the 
dispersive regime. In Table HT| we list the gate times, the 
decoherence rates, and the ratios of gate times to deco- 
herence times for the gates discussed above at k = 4 MHz. 



One can estimate the gate fidelity approximately as 

F = e- T °/ T « w 1 - T g /r d (24) 

which depends on the ratio between the gate time and the 
decoherence time. Using the parameters given above, we 
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k (MHz) 

Figure 3: Gate fidelity versus damping rate K by numeri- 
cal simulation. Solid curve: SWAP gate and dash-dot curve: 
Hadamard gate. Two fluctuators are included in the simu- 
lation for both gates. For the Hadamard gate, we use the 
parameters A 2 = -2tt x 80 MHz and g 2 — 2ir x 30 MHz for 
TLS n = 2. 



can estimate the ratio Tg/r^. In Table ITTl it is shown that 
T g/ T d ~ 10~ 2 for the two-qubit gate in our scheme and 
T g/ T d ~ 2 x 10 -2 for the SWAP operation in the Cirac- 
Zoller gate at the damping rate k = 4 MHz. The gate 
fidelities can thus reach 0.99 with our protocols. This in- 
dicates that the Josephson junction resonator is an effec- 
tive tool in generating high fidelity quantum operations 
even at a quality factor (Q = u> c / «) of only Q — 7800 and 
hence can be used to demonstrate quantum coherence in 
the TLS fluctuators. 

The term a nz (a + aJ) in the residual coupling H x in 
Eq. (frTj) can induce additional decoherence. Here, quan- 
tum fluctuations of the resonator mode cause extra noise 
even at zero resonator amplitude. The spectrum of the 
quantum fluctuations can be derived as 



(aa% = ((oj - lo c ) 2 + k 2 /4)" 



from which the decoherence rate can be derived as 



(25) 



(26) 



As can be seen, the decoherence rate is to the 4th order 
of g n /A nc when the driving amplitude e is comparable 
with the detunings. With the above parameters, we find 
that t^ 1 sa 1 kHz which can be neglected during the gate 
operations. 

To study the effectiveness of the quantum logic gates, 
we perform numerical simulations on the gate opera- 
tions using the full Hamiltonian H t . The decay of the 
resonator is simulated using the Lindblad master equa- 
tion 2 ^. This simulation includes both the effect of the 
residual coupling in Eq. (fl7f and the effect of the res- 
onator decoherence. We calculate the fidelities of the 
Hadamard and SWAP gates using the method prescribed 
by Nielsen's formula for the gate fidelity 2 ^, over a wide 
range of resonator decay rate. The results are plotted 



in Fig. [3l This simulation shows that the fidelity can be 
higher than F — 0.99 for n < 5 MHz for single-qubit and 
two-qubits gates, which also agrees with our estimations 
above. Hence, at moderate decay rate for the Josephson 
junction resonator, high fidelity quantum logic gates can 
be achieved. 



VI. DISCUSSIONS 

As demonstrated in recent works, the superconducting 
resonator can be used to detect the states of qubits or 
fluctuators^^. In the strong damping regime, the ampli- 
tude of the superconducting resonator adiabatically fol- 
lows the dynamics of the TLS fluctuators, as was shown 
in our previous work* 7 . A phase sensitive measurement of 
the resonator can hence provide a direct measurement of 
the TLS states. In the moderate damping regime where 
the damping rate is weaker than the coupling constant, 
a measurement of a TLS can be performed by adjusting 
the resonator frequency to be in the vicinity of this TLS, 
but with the condition g n <C |A nc |. Here, the Stark shift 
resulting from this TLS is much stronger than that from 
other qubits. A measurement of the resonances in the 
transmission spectrum of the resonator can then provide 
a measurement of the qubit stated. Such measurements 
can be realized with current electronics where the res- 
onator and the driving can be adjusted and switched on- 
and-off in nanoseconds, much faster than the decoherence 
time of the fluctuators. 

Our scheme can also be extended to TLS fluctuators 
inside different junctions. Because the wavelength of 
the microwave mode of the superconducting resonator 
is much longer than the dimension of this circuit, fluctu- 
ators in different junctions can couple to the same res- 
onator mode when the junctions are connected by a su- 
perconducting loop. As is illustrated in Fig. [T](b), two 
junctions are connected to the central superconducting 
island labelled by $ that is associated with the junc- 
tion resonator. It can be shown that effective coupling 
between fluctuators in the two junctions can be derived 
exactly as described by Eq. (fTTjl . In this configuration, 
quantum logic gates can be performed with essentially 
the same approach as was presented above. This circuit 
can also be extended to include multiple junctions. This 
system is thus intrinsically "scalable" where fluctuators in 
multiple junctions couple nonlocally. Note that the fre- 
quency of the superconducting resonator is determined 
by the total capacitance, the total effective Josephson 
energy, and the inductance in the circuit, and will de- 
crease as the number of junctions increases. This can set 
a limit on the number of junctions that can be connected 
into the circuit. 

TLS fluctuators have been studied for a long time. Pre- 
viously, the fluctuators are often considered as a source 
of decoherence in superconducting qubits, causing the so- 
called 1// noise. In this work, we focused on studying 
the coherent manipulation of the fluctuators which can 
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provide insights about the dynamics, the coupling mech- 
anism, and the relaxation of the fluctuators in supercon- 
ducting devices. Although the success in implementing 
universal quantum logic gates makes the TLS fluctuators 
potential candidate for quantum computing, we want to 
emphasize that the main aim of this work is to provide a 
practical scheme to demonstrate the coherence behavior 
of the fluctuators. Our scheme can be useful for current 
experiments that investigate the coupling between the 
fluctuators and superconducting resonators^ ^ 14 ' 17 ' 18 . 

VII. CONCLUSIONS 

To conclude, we have shown that universal quantum 
logic gates can be implemented on spurious TLS fluctu- 



ators via the coupling between the fluctuators and the 
Josephson junction resonator. Taking into account the 
full Hamiltonian of the coupled system and the effect 
of the noise, our numerical simulation of the quantum 
operations showed that quantum logic gates can be per- 
formed with high fidelity even at resonator decay rates of 
a few megahertz. We have used practical parameters for 
the junction resonators and the fluctuators in this study. 
Our work hence indicates that quantum coherence and 
quantum manipulation of TLS fluctuators can be read- 
ily demonstrated. The results here can be generalized 
to other types of superconducting resonators which are 
explored in recent experiments^ 4 -. 



* Electronic address: ltian@ucmerced.edu 
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